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ABSTRACT 
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Degree and Year: Master of Science, 1989 

Thesis directed by: Dr. Kevin J. Webb 

Assistant Professor 

Electrical Engineering Department 

A frequency selective surface (FSS) composed of apertures in a metallic 
sheet is known as an inductive FSS. The infinite inductive FSS theory is de- 
rived and the aperture fields are solved by a spectral domain formulation with 
method of moments solution. Both full domain and subsectional basis func- 
tions are studied. A locally planar technique (LPT) is used to determine the 
forward scattered field from a generally shaped inductive FSS with arbitrary 
illumination. 

An experimental FSS study describing test methods and results is pre- 
sented to verify both the infinite FSS theory and determine the accuracy of the 
LPT. The infinite planar FSS transmission coefficient measured as a function 
of frequency is compared with theory. The inductive FSS s forward scattered 
field is experimentally measured and compared to the LPT. A series of experi- 
ments with the spherical feed source at varying distances from the planar FSS 
is used to determine the limitations of the LPT in terms of the incident phase 
error over a unit cell. The limitations of the LPT due to surface curvature 
is investigated by an experimental study of the scattered fields from a set of 


preceding page blank not filmed 


hyperbolic cylinders of different curvatures. Prom the experiments, guidelines 
for applying the LPT are developed. 


iv 


TABLE OF CONTENTS 


CHAPTER 1. INTRODUCTION 

CHAPTER 2. INFINITE INDUCTIVE FSS ANALYSIS 

2.1 Introduction 

2.2 Spectral Domain Formulation 

2.3 Spectral Moment Method 

2.4 Basis Functions 

2.5 Transmission Coefficient 


CHAPTER 3. LOCALLY PLANAR TECHNIQUE (LPT) 

3.1 Introduction 

3.2 Region I. Feed Source 

3.3 Region II. Inductive FSS 

3.4 Region III. Forward Scattered Fields .... 

3.5 Locally Planar Technique Summary .... 


CHAPTER 4. FREQUENCY RESPONSE OF THE INDUCTIVE FSS: 
NUMERICAL AND EXPERIMENTAL RESULTS 

4.1 Inductive FSS Design 

4.2 Transmission Coefficient Magnitude: 

Numerical and Experimental Results 


V 


4.3 Transmission Coefficient Phase: 


Numerical and Experimental Results 49 

CHAPTER 5. INDUCTIVE FSS FORWARD SCATTERED FIELD 
EXPERIMENTS AND MODELING 

5.1 Introduction 53 

5.2 Scattered Field Experiments With Plane Wave Incidence 

on a Planar FSS 56 

5.3 Study of the LPT Accuracy as a Function of Incident 

Field Curvature 65 

5.4 Study of the LPT Accuracy as a Function of Surface 

Curvature 78 

CHAPTER 6. CONCLUSIONS 86 

REFERENCES 89 


vi 


/ 


ACKNOWLEDGEMENTS 


The author would like to thank his thesis advisor, Dr. Kevin J . Webb, 
for his guidance and advice in the FSS study. A very special thanks goes 
to Billy B. Williams of the RF Systems Section for the countless hours of 
effort in performing the numerous FSS experiments. The author also thanks 
Preston W. Grounds and Richard F. Schmidt for their helpful explanations 
and clarifications of the FSS theory. 








CHAPTER 1 


INTRODUCTION 

The frequency selective surface (FSS) is a useful filter for electromagnetic 
radiation ranging from microwave to optical frequencies [1]. Two types of 
FSSs exist. The capacitive type, illustrated in Fig. 1.1, is composed of metal- 
lic patches on a dielectric substrate. The inductive type, shown in Fig. 1.2, 
consists of apertures in a metallic sheet. The dots in both figures signify that 
the patches and apertures are a part of an infinite FSS. The response of the 
infinite planar inductive FSS is determined from a spectral moment method 
solution derived in Chapter 2. A typical capacitive and inductive FSS response 
about the resonant frequency, /o, is shown on the right side of Figs. 1.1 and 
1.2. Note the magnitude of the reflection coefficient in the capacitive case is 
equivalent to the magnitude of the transmission coefficient for the inductive 
case. This is valid if the patches rind apertures can be interchanged to trans- 
form the one type of FSS to the other. The duality of the frequency responses 
can be explained by Babinet’s principle [2j. 

In many FSS applications, such as radomes and subreflectors, the FSS 
is curved and has non-planar illumination. Figure 1.3 is an example of a 
reflector antenna utilizing a hyperbolic FSS to separate the S-band Cassegrain 
feed and K-band prime feed [3]. The FSS shape is a hyperboloid of revolution 
and the incident wave is spherical. The general analysis of curved FSSs is 
an extremely difficult problem. An approximate method, called the locally 
planar technique, (LPT) can be used to predict the response from a planar or 
curved FSS with arbitrary illumination. This technique has been applied to 
the case of a parabolic cylinder FSS by Ko and Mittra [4]. This technique, 
presented in Chapter 3, breaks the FSS into a number of subarrays each of 
which is assumed to be a segment of an infinite planar surface. The infinite 
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Figure 1.1. Capacitive FSS metallic patch geometry and the magnitude of the 
reflection coefficient as a function of frequency. 



Figure 1.2. Inductive FSS apertures in a metallic sheet and the magnitude of 
the transmission coefficient as a function of frequency. 
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Figure 1.3. Dual reflector antenna system with an inductive FSS subreflector 
separating the prime and Cassegrain feeds. 
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FSS theory is applied to each subarray and the total response is found from a 
superposition of individual responses. 

However, there are limitations for the FSS geometries which can be solved 
by the LPT. At some value of incidence phase front curvature and or surface 
curvature, the LPT is no longer valid. In the literature, the limitations of the 
LPT for predicting the scattered fields from the general FSS problem has not 
been addressed. Thus, a major goal of the FSS study was to determine the 
regions of validity for the LPT theory. To study the LPT as well as verify 
the infinite planar FSS analysis, an inductive FSS was designed and fabri- 
cated. An extensive experimental FSS study was undertaken at the NASA / 
Goddard Space Flight Center. Chapter 4 describes the infinite planar experi- 
mental techniques and results from the transmission coefficient measurements 
as a function of frequency. In Chapter 5 two types of LPT experiments and 
scattered field pattern results are presented. The first LPT experiment is for 
a planar FSS with varying incident field curvature. The second LPT experi- 
ment is for a set of hyperbolic cylinders of varying curvature with near planar 
illumination. By comparing theory and experiment, guidelines for applying 
the LPT are developed. From the guidelines, a given FSS problem can be 
evaluated for its suitability for LPT analysis. 
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CHAPTER 2 


INFINITE INDUCTIVE FSS ANALYSIS 

2.1 Introduction 

The infinite inductive and capacitive FSS has been analyzed by a number 
of authors over the years [4]- [8]. The spectral domain moment method tech- 
nique applied to the inductive FSS in this chapter is based primarily on the 
work of Tsao and Mittra [8]. A good summary of the technique is given in 
the text titled The Spectral Domain in Electromagnetics [9]. The FSS theory 
presented here is for an infinitesimally thin FSS illuminated by a plane wave. 
The analysis begins from basic electromagnetic theory. From a spectral domain 
method of moments solution, the unknown aperture fields are calculated. 

2.2 Spectral Domain Formulation 

The analysis begins by considering the inductive FSS illustrated in 
Fig. 2.1. The FSS is composed of rectangular apertures placed in a rectangu- 
lar grid arrangement. A unit cell of the FSS is defined as the region enclosed 
by the dashed line in Fig. 2.1. The dimensions of the unit cell are equal to 
the aperture-to-aperture spacing in each direction. Figure 2.2 illustrates the 
incident plane wave geometry. The incident electric and magnetic fields are 
given by E ,nc and H‘ nc . The plane wave’s propagation vector, k ,nc , is oriented 
at an angle 6 from the z-axis and at an angle (f> from the i-axis. The analysis 
assumes the FSS to have zero thickness. The experimental FSS had a thickness 
much less than the smallest wavelength of interest, so the assumption of zero 
thickness is valid in the theory. The scattered magnetic field from an aperture, 
H 4 , is derived from the electric vector potential F 

H* = -jue F + -r~V(V • F). (2.1) 

JUfi 
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Figure 2.1. Inductive FSS aperture geometry. The FSS shown is a small portion 
of the larger “infinite” FSS. 
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Figure 2.2. Plane wave incidence geometry on the inductive FSS. 
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The electric vector potential, F, is found from the magnetic current K. K 
in turn is derived from the equivalence principle for electric fields [10]. The 
equivalent magnetic current is placed above a perfect electric conductor (PEC) 
and has the value 


K = E x n. (2.2) 

The orientation of the unit vector n is normal to the FSS in the direction of 
the half space in which the fields are desired. The electric vector potential F 
in the (x, y) plane, based on the definition E = —V x F, is derived from K 
and the Green’s function G by 


F(x,y) 


(F x \ _ f G * K x \ 


(2.3) 


The symbol * is the convolution operator. K x and K v are the components of 
K in the a x and a y directions respectively. The Green’s function G is defined 
by 


G(x,y) 


e~ jk ° r 

4nr 


(2.4) 


In the (x,y) plane, 


r = \/x 2 + y 2 . 


(2.5) 


The free space propagation constant is 


2n 2irf 

ko = “T" = • 

A c 


( 2 . 6 ) 


The frequency of operation is /, and c is the speed of light. Using the identity 
dyadic /, G is defined as 
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(2.7) 


g = g/ = g(£ ;). 

Equation (2.3), for the electric vector potential, can be written as 

F(x,y) = G*K. (2.8) 

Returning to (2.1) for the scattered magnetic field and expanding in the 
Cartesian coordinate frame yields 

H*(x,y) = —ju>e(F x a x + F v a y ) + 


1 

jufi 


d i 

(dF x 

dF' 

, d / 

<dF x 

8F t \ 


\ dx 

4- y 

dy j 

) ai + ^( 

^ dx 

+ ~aF/ \ 


(2.9) 


Equation (2.9) can be written into matrix format and simplified by pulling out 
the term and noting = ^^e scattered magnetic field in the 

(x, y ) plane is 


i /a* . jfc2 A. A. \ 

H'(x,y) = - — ( Sx 8 e ° a’’ ’.2 ) F (*'V)- 

)^\ &■§; fc + k SJ 


( 2 . 10 ) 


H* is transformed into the spectral domain by the Fourier transform defined 
by 


CO OO 

H '(<*,/?)= J /H*(x,y)e-™dxdy, (2.11) 

— oo — OO 

where H J denotes the transform of H*. The spectral variables a and (3 take 
on discrete values, the phase constant values for the Flouquet modes, which 
for the rectangular cell Eire [11] 
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+ K nc 


771 = 0, ±1, ±2, 


(2.12a) 


2irm 


— 


/ h = + ^ nc n = 0,±l,±2,--. (2.126) 

c v 

tj." c and t* nc axe the projections of the incident propagation vector into the 
(x, y) plane. The constants a and c are the unit cell dimensions in the x and 
y directions respectively. From Fig. 2.2, 


tj. nc = to sin 6 cos (f> 


(2.13a) 


k' y nc = to sin 6 sin <f>. 

(2.136) 

The equations derived in the spatial domain are transformed into the spectral 
domain by means of the following equivalent operators: 

A(x, y) * B(x, y) =$> A(a, p)B(a, P) 

(2.14a) 

— A(x, y) =» ja A(a, p) 

(2.146) 

A(x, y) j(3 A (a, p). 

oy 

(2.14c) 

The transform of the scattered magnetic field, H 4 , is found by a 
and (2.14c) to (2.10). The result is 

pplying (2.14b) 

*•«*•» -sc? C-'«; 

(2.15) 
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Applying relation (2.14a) to the electric vector potential equation (2.8) results 
in 

F (aj) = GK (2.16) 

The Fourier transform of the Green’s function, is [12] 


G(aJ) = 




2 \Ao ~ a 2 - fi 2 

for the propagating modes, and 


* 2 >(« 2 +/? 2 ) 


(2.17a) 


G(a,(3) = 


k 2 0 < (« 2 + P) 


2^a 2 TW~ K 

for the evanescent modes. Substituting (2.16) into (2.15) yields 


(2.176) 


H* 




(2.18) 


The scattered magnetic field in the spatial domain, H a , is found by taking the 
inverse Fourier transform, giving 


oo oo 

H* = J J H*(o, f3)e 3 ( ax+ P y ^ dad/?. (2.19) 


— oo — OO 


The integral becomes a summation since the magnetic current, K, is non-zero 
at discrete points only. H J in the spatial domain is 


H' 


1 ^ — a 2 —a/3 ^ (^i^ e j(o,z+Py) 

L 2. ^ -a/? ^o 2 -/?V 

J ' m=— oon=— oo 


( 2 . 20 ) 


To uniquely determine the magnetic currents in the aperture, the bound- 
ary conditions at the z = 0 plane will be evaluated. Figure 2.3(a) illustrates 
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(b) Equivalent magnetic field representation. 

(c) Method of images fields and currents in free space. 



the original problem with the aperture field E and the incident magnetic 
field in the plane of the FSS, Hj. n „ c . Figure 2.3(b). illustrates the equivalent 
magnetic current backed by a PEC for both incident and transmission regions. 
In the incident region, the equivalent magnetic current K and magnetic field 
H ,nc are located an infinitesimal distance e above the PEC. From the method 

xy 

of images, when the PEC is removed, equivalent fields and currents can be 
determined for the free space region illustrated in Fig. 2.3(c). The resulting 
magnetic field and magnetic current have doubled. 

The magnetic field in the (x,y) plane due to the magnetic current K is 
given by H*. The total magnetic field, H + , in the incident region at z = +e, 
is 


H + = 2H£ c + 2H\ (2.21) 

In the transmission region, for the continuity of the tangential electric field, 
the magnetic current K — = — K~^. The magnetic field in the transmission 
region, at z = — f , shown in Fig. 2.2c, is due entirely to 2K and given by 


H~ = -2H'. 

The aperture field is continuous, so as e — * 0 


( 2 . 22 ) 


H" = H + . (2-23) 

Equating (2.21) and (2.23) and simplifying yields the relation between the 
incident and scattered field on the (x, y ) plane 


H a = 


Hj w w c (x,y) 

2 


(2.24) 


Substituting the above into (2.20), replacing K + with Exa„ and simplifying 
yields 
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H 


»r»c 

xy 


i 

CJ/i27T 2 


OO OO s 

E E ( 

m— -oon=- oo ' 


a/? 

P 2 ~ K 


k ° GEe j(ax+fly \ 

-**P ) 


(2.25) 


the inductive FSS equation. In the next section, a solution technique will be 
presented to solve for the aperture electric fields. 

2.3 Spectral Moment Method 

In this section, the basis function weights will be found by the moment 
method implemented in the spectral domain. The inductive FSS equation 
(2.25), relating the unknown aperture field to the known incident field has the 
form 

L(f) = g. (2.26) 


L is given by 


L = 


. 7 


OO OO s 

E E ( 

m= — oon= — oo ' 


&P 

P-kl 


*1 ^ Q e j(<*z+0v) 

-<*P J 


(2.27) 


The unknown f is defined as 


f =(!:), (2.28) 

where 

N, 

E z = (2.29a) 

n= 1 

E, = (2.294) 

n=l 

E Xn and E Vn are the basis functions described in Section 2.4. 
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In (2.26) g is given by 


g 



(2.30) 


The unknown complex constants, c„, are determined by constructing a 
set of N = N x + Ny linearly independent equations. The equation set is formed 
by an inner product defined by 


< A - B > = // 


(A x B*) • a* ds, 

aper 


(2.31) 


where the symbol * denotes the complex conjugate. The inner product is 
taken on both sides of equation (2.26) with respect to the testing functions b n 


f;<L(f),b n ) = (g,b„). (2.32) 

n=l 

In the Galerkin’s method, the testing functions are chosen to equal the basis 
functions [13], so 


b n 6 {£*„a*,£ y „a*} ( 2 - 33 ) 

Equation (2.32) is composed of an x-component and y- component given by: 

OO CO 

3 - Y Y G(aPE x + (kl-a 2 )E y )e* ax+ ™ = tf* nc , (2.34a) 

uu,2n 2 

' m= — oon=- oo 

co co 

— I — - y y G({f -kl)E x -aPEy) e* ax+M = H { y nc , (2Mb) 

where E x and E y are the Fourier transforms of the aperture electric field in 
the a x and a y directions, respectively. To form a compact set of equations, 
the following substitutions will be made: 
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X (a/?) 

(2.35a) 


(2.356) 


(2.35c) 


The N simultaneous equations axe formed by computing the inner products 
between (2.34) and the testing functions. Inner products are taken between 
the E x testing functions and (2.34b) and between the E y testing functions and 
(2.34a). This is a direct consequence of the cross product in the inner product 
definition. The form of the inner product on the left-hand side, (LHS} x , using 
(2.34b) for the x-directed testing function is 


(LHS), 


, oo oo 

= ( E E 

' m= — oo n— — c* 


( ^3 E x q\Ey)e 


j(ax+0y) 


» 


E x 


-> 


(2.36a) 


oo oo 

= (a, X a.) ■ a, / / £ £ (* A - q,E,)El n dxdy. 

J J opcr m= _ 00 n=; _ 00 

(2.366) 

The integral can be written as a Fourier transform, of E* n . The inner product 
in (2.36) becomes 


OO OO 

(MS), = - £ £ (93 E,-q,E,)E;„. 

m = — oo n= — oo 


A similar result for the E y testing functions is 

oo oo 

(LHS) y = Y, E (toE* + toE 9 )E; n . 

m= — oo n— — oo 


(2.37a) 


(2.376) 
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The y-component of the incident magnetic field on the (x, y ) plane is given by 

jjinc _ ^•inc e >(or(0)x + /?(0)») (2.38) 

The inner product on the right side for the x-directed testing function can be 
written as the zeroth spectral component of a Fourier transform 

(MS), = <Hjr,E,„> = -W"‘E;Ja(0),m). (2.39) 

The similar result for the y-directed testing functions is 

(RHS) y = H[ nc E* yn {a{ 0),/J(0)). (2.40) 

The final set of equations for the inductive FSS is 


oo oo 

E E 

m=— oo n= — oo 


oo oo 

y y + qiE y )E* VN f 

m=— oo n= — oo 

E E (»& - 

m = — oo n= — oo 


oo oo 

y y (q^Ez — q\E y )E XN _ 

m= — oo n= — oo 

Substituting in the basis function expansion for the electric field (2.29) and 
writing in matrix form yields an equation of the form 

[A][C] = [ft]. (2.42) 








(2.41) 
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The matrix [A] is given in Fig. 2.4. For computer implementation, the inner 
product summations are taken over a finite m and n set. The number of 
basis functions and spectral terms for a satisfactory convergent solution will 
be discussed in Chapter 4. The vector [B] is 


(“(»)«»)) 


[fl] = 


(n(0)/J(0)) 


\b" c e;„Jo{ o)m)J 

Vector [ C ] is the unknown weighting constants for the aperture field 


(2.43) 


[C] = 


( \ 


C V\ 


\ C VN, / 


(2.44) 


The unknown [C] in (2.42) is solved for by a numerical technique which is 
equivalent to 


[C] = [A]-’[B]. (2.45) 

... / 

The aperture fields are found by a superposition of the weighted basis func- 
tions defined in equation (2.29). The aperture field will be used in Section 
2.5 to determine the amplitude and phase of the transmission coefficient. In 
Chapter 3, the forward scattered field from a finite FSS will be found from the 
transforms of the aperture fields. 
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Figure 2.4. Matrix [A]. 


2.4 Basis Functions 


Two types of basis functions, the full domain and subsectional, are stud- 
ied. The desired basis function set is one which can model the aperture fields 
with the least number of unknown weights, c n . The basis functions ideally 
suited for the problem of basic geometrical aperture shapes are the full do- 
main basis functions where apriori knowledge of the edge condition can be 
included. The experimental FSS described in Chapter 4 has rectangular aper- 
tures ideally suited for this type of basis function. The full domain basis 
functions have the property that each basis function spans the entire aperture. 
The subsectional basis functions, on the other hand, define each basis function 
over a different area of the FSS. The subsectional basis functions axe useful for 
modeling the aperture fields for general aperture shapes. However, the sub- 
sectional method requires a large number of basis functions to model the fields 
accurately. Hence, the CPU time is much greater than for the full domain 
approach. In this section, both basis sets will be defined over a rectangular 
aperture and the Fourier transforms will be calculated. The basis functions 
are derived for the x-directed electric field only since the y-directed basis set 
is similar. 

Full Domain Basis Functions 

The basis functions used to model the x-component of the aperture electric 
fields in a rectangular aperture of dimension 6 x d are [14] 


where 



(2.47a) 


p = 0, 1,2, and q = 1,2, — 


(2.476) 


Tp(x) is the Chebyshev function of order p. The first few Chebyshev functions 
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are given by [15] 


T 0 (x) = 1 
Ti(x) = x 
T 2 (x) = 2x 2 - 1. 


(2.48) 


Figure 2.5 shows the shape of the major full domain basis functions. The 
result of the Meixner edge condition is seen at x = ±|. 

The Fourier transform, (2.11), of the basis function is found by separating 
the functions into x and y components and transforming each into the spectral 
domain. In the space domain let 


E In — f x (*)/y (y) a i- 


(2.49) 


The function f x is defined by 

T (iz) 

/,(*) = P,(x) V (2-50) 

(>-(*)) 


p b (x) is the unit pulse of width b centered at the origin. The transform of f x 
is 


/,= /p l( x)T p ( f)(l-(f) 2 ) ( 251 ) 

— oo 

Msiking the change of variables x\ = ^ and noting the integral is limited by 
the pulse yields 


l 

/, = r,(x!)(l - (xj) 2 )“ i e-- ,( "* I,/2> dxi. (2.52) 

-1 
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Figure 2.5. Full domain basis functions for modeling the E* aperture electric 
field. 
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The integral is solved by associating terms with the known integral [16] 


J T„(i)(l - i 2 r*e->"' dt = (-l)"0T<rJ«M, (2.53) 

-1 

where J n (w) is the Bessel function of the first kind of order n. The transform 
/ z becomes 

h ( 2 - 54 ) 

The y-dependent term is 

fy = Pd(y) sin (y + 0 ^ • ( 2 - 55 ) 

The transform is taken by separately transforming the sine and pulse func- 
tions and convolving the two functions in the spectral domain. The Fourier 
transforms needed are 


Pd(y ) =S>dsinc^y^ 

in (t ( y + 0) ^ “ 3 * - T ~> " + t)] e ’ M2 

oo 

h{y)My) =► ^ J fi(y)MP - 0 d *- 

— oo 


(2.56) 


The 8 function is defined by 

*>-{r (2 - 5?) 

The convolution integral is simplified by the 8 function. The resulting expres- 
sion for f y is 
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To simplify (2.60), note that the exponentials can be represented by 

e”*/ 2 = (j) q 
e -jg*/2 = 


(2.58) 


(2.59) 


The expression for f y becomes 


f -JA 

Jy n 


c j y sinc (Jj- - y) ~ sinc (^r + ?) • ( 2 - 6 °) 


The two transform equations (2.54) and (2.60) are combined to yield E Xn : 

e„ =^or +1 (-i)'j,(f ) 

>-)•-“ (t - t) -(-•»■)• -“(t + t)]- 

Similarly the result for E Vn is 

or sinc ( y - q -f) - (-j)' sinc (t+t)]- 


Subsectional Basis Functions 

The subsectional basis functions used are called rooftop functions since 
they are composed of triangles in one direction and pulses in the orthogonal 
direction [17], [18]. As was the case with the full domain basis functions, only 
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the E* basis function will be discussed. Figure 2.6(a) shows a single rooftop 
basis function located at a position (x„,y n ) in the aperture. As was dis- 
cussed in the full sectional case, the electric field in the x-direction goes to 0 
at y = ±| and peaks at x = f . Thus the rooftop functions for E x are placed 
with the roof line parallel to the x axis to best model the field. Figure 2.6(b) 
shows the layout of the pulses spaced D x apart. N x is the total number of 
pulse functions in the x-direction. Figure 2.6(c) illustrates the N y overlapping 
triangles each having a base width of D y . The pulse and triangle functions 
P D io( z ) and &D,y 0 (y) centers x 0 and y 0 are defined as 


Pd,x o 



0 , 

1 , 


for |x — Xo| > y', 
for |x — xo| < y- 


&D,y o(y) 



for |y — yo 1 > y; 
for | y - y 0 | < y ■ 


A single rooftop electric field basis function is defined by 


(2.62) 


(2.63) 


E Xn (x,y) = PD',r„(x)A D ,,y n (y)- (2.64) 


The pulse and triangle base widths are 


D x = 


b 

N~ z 


and 


D y 


2d 

N y + 1' 


(2.65) 


The basis number, n, is related to the rooftop basis function center (x n ,y„) 
by a basis function layout scheme. 

The Fourier transform of the basis function (2.64) is obtained from the 
transforms of the centered pulse and triangle combined with the shift property 
in the space domain. The transform relations of interest are 
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(b) (c) 


Figure 2.6. Subsectional basis functions for the E x aperture field. 

(a) Single rooftop function with center (x n ,j/ n ). 

(b) Pulse function layout. 

(c) Overlapping triangle functions. 
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( 2 . 66 ) 


faD x \ 

Pd., o(x) => D t sine I I 

a -Py • 2 / PDy \ 

A D , )0 =»-^smc J 

/(* -*o) = /(a)e»'“- 

/(»-»>) = /w^ w - 

From the relations given in (2.66), the Fourier transform of E Xn is 




D x D y . / al? 

— - — sine — — 

2 V 2 




,>(ar n +0y n ) 


(2.67) 


The above equation can be used to find the Ey basis functions if the following 
substitutions axe made: x y, a b *-+ d, and n x *-* n y . 

2.5 Transmission Coefficient 

The planar FSS’s transmission coefficient as a function of frequency and 
incident angle is a principal means of characterizing the behavior of the FSS. 
In typical operating modes, a single plane wave is transmitted. However, 
at some frequency above resonance, and or at large incident angles 0, other 
transmitted modes appear. This effect, which can be compared to grating 
lobes in a phased array, is called the Wood’s Anomaly. Thus, in general a 
number of transmission coefficients, corresponding to each propagating mode, 
must be calculated. The magnitude of the transmission coefficient will be 
found by first defining the total transmission coefficient in the spatial domain. 
Then, by invoking Parseval’s theorem the power transmission coefficient of 
each propagating wave will be determined. 

Referring to Fig. 2.2, the magnitude of the total transmission coefficient, 
|T|, is defined as the power transmitted normally through the FSS, Ptrans, 
divided by the power incident normally on the surface, Pine- The expression 
for |T| is 
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( 2 . 68 ) 


Re// apgr (E a xH a )-(-a z )da 
11 Re / /„„« (E ,nc x H*" c ) • (— a z )da" 

cell 

E a and H a are the toted fields in the aperture related by 

— jufiHa = V x E a . (2.69) 

From the above equation, the magnitudes of the aperture fields are related by 
the impedance of free space, Zo 


|E 0 |= Z 0 |H a |. (2.70) 

The numerator of (2.68) can be simplified by (2.70) to yield 

A £ 

Ptrans = / / 4" I E a | 2 COS 9 dx dy. (2.71) 

J=± J-i_ z 0 

The total aperture electric field, E a , is composed of the electric field in the 
plane of the aperture, E, and an unknown E z component. Given that the 
propagation vector for the radiated plane wave is k p , the electric field must 
satisfy 

E a • k p = 0 (2.72) 


E z can be found as 


E z 


-(E x kP + E y kP) 


(2.73) 


The integral (2.71) can be simplified by the Parseval's relation for periodic 
signals [19] 
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(2.74) 


XoY 0 


f j I A(x, y) dxdy = ZEK«))f, 

J Xn JYo m n 


where Xq and Yq correspond to the period of the function in the x and y- 
directions, respectively. The integral in (2.71) can be simplified by Parseval’s 
theorem with Xq = o. and Yq — c. The result is 


Ptran» = — AC COS 
6 0 


i o ^2 |E a ^a(m p ),^(n p ))| , (2.75) 


( m p » n J») 


where 6 is the angle between the propagation vector of a given plane wave and 
the -z-axisj and the summation is taken for pairs of ( rripyUp ) corresponding to 
propagating modes. The propagating modes are found from the field variation 
along the z-axis, e~i k,z . The propagation constant k z must be positive for 
power to be transmitted, where k z is defined as 


k z = \fk 2 — a 2 (m) — (2.76) 

Each propagating mode must have the pair ( m p ,n p ) satisfying the inequality 


k 2 > a 2 (m p ) + /? 2 (n p ). (2.77) 

Returning to (2.68), the denominator, P inc , represents the power incident on 
a unit cell. From the incident geometry, shown in Fig. 2.2, 


P inc = ac | E ,nc 1 1 H ,nc | cos 9 = Z 0 ac | II inc | 2 cos 0. (2.78) 


The power transmission coefficient for each transmitted mode found by 
dividing (2.75) by (2.78) is 
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( 2 . 79 ) 


I T(m p ,n p ) 


a (a(m p ),^(n p ))| 
Zl | H ,nc | 2 


The phase of the transmission coefficient is the phase change of the inci- 
dent wave upon transmission through the FSS in a direction given by a(m p ) 
and 0(n p ). If $ is the phase of the incident plane wave and ^ is the phase of 
the Fourier transform of the aperture electric field in the direction of interest, 
then the phase of the transmission coefficient is $ + 
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CHAPTER 3 


LOCALLY PLANAR TECHNIQUE 
3.1 Introduction 

The LPT, introduced in Chapter 1, is an approximate technique which 
applies the infinite FSS theory to general FSS shapes with arbitrary incidence 
subject to some constraints. Before beginning the LPT study, it is worth 
mentioning the types of FSS problems which have been numerically solved to 
date. The finite FSS problem can be solved by calculating each aperture’s 
field distribution by matrix inversion for a small number of apertures or by 
iterative techniques for a large number of apertures [20]. The analysis of the 
general curved FSS is quite difficult. The curved FSS analysis presented in 
the literature has been for the case of patches located on a cylindrical surface 
[21], [22], The numerical analysis for problems like the hyperbolic subreflector 
is extremely complex and it appears the LPT may provide the only feasible 
solution for this type of FSS problem. 

In this chapter, the LPT will be developed by applying the theory to the 
problem of determining the forward scattered field pattern from a generally 
shaped inductive FSS with arbitrary illumination. Figure 3.1 illustrates the 
overall geometry and defines the coordinates. In Region I, the incident fields 
on the FSS are determined from a feed source model. In Region II, the FSS is 
broken up into N subarrays, each of which is considered to be locally planar. 
The coordinate system (X’(n), Y^n), Z’(n)) is the local coordinate frame for 
each subarray specified by n. Invoking the infinite FSS theory at each subarray 
yields the approximate aperture electric field. The forward scattered far fields 
are determined in Region III from the Fourier transforms of the aperture fields. 
The scattered field at the point (x s ,y s , z s ) is found by summing up the fields 
radiated from each subarray. In Chapter 5, the LPT will be compared with 
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Figure 3.1. Locally planar technique geometry overview for the n th subarray. 
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experimental results to determine the limitations of the theory. 

3.2. Region I. - Feed Source 

The feed source radiates a linearly polarized field with the electric 
field oriented parallel to the x-axis. The feed phase center is located at 
(x For the experiment, the feed chosen was a pyramidal horn with 
equal beamwidths in the x-z and y-z planes. The experimental feed pattern 
was closely modeled in amplitude by a c,os m 9\ cos m 62 approximation. Fig- 
ure 3.2 illustrates the feed geometry. The calculated and measured feed pat- 
tern for a value of m = 2.355 is shown in Fig. 3.3. The solid line is the cosine 
numerical approximation. The dashed line and dotted line are the H-plane and 
E-plane measured patterns, respectively. In all experiments, the FSS is in the 
far field of the transmit horn. Thus, the magnitude and phase of the incident 
field at each subarray n on the FSS can be modeled by 

E inc (n) = — r cos m 0\ cos m d 2 e~ ikoR{n \ (3.1) 

J?(n) 

R(n) is the distance from the feed phase center to the location of each subarray. 
The incident electric field in the global coordinate frame, (X, V, Z'), is E (n) 
where 

E x = E inc (n ) cos l?i 
E y = E' nc (n ) sin 0i sin 02 
E z = — E ,nc (n) sin 0i cos0 2 . 

The components of the incident propagation vector are 

k x = k 0 sin 0i 
k y = k 0 cos 0i sin 0 2 
k z = ko cos 0i cos 0 2 . 


(3.2) 


(3.3) 
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Figure 3.2. Pyramidal horn feed geometry. 
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Figure 3.3. Pyramidal horn approximate numerical and experimental patterns. 


34 



From Fig. 3.2 the incident magnetic field is 


H ,nc (n) 


k ,nc (n) x E ,nc (n) 
Zq ko 


(3.4) 


where Zq is the impedance of free space. 

3.3. Region II. - Inductive FSS 

The general shaped FSS is first divided into N subarrays. Each subarray 
is made up of one or more apertures. Ideally, the subarray should be a single 
aperture in size. This is because the smaller the subarray, the more accurately 
the surface geometry and incident field contour is represented. The main 
reason to use subarrays composed of many apertures is to reduce the number 
of infinite FSS calculations, and hence the CPU time. In the analysis which 
follows, the subarray is composed of a single aperture and N represents the 
number of apertures in the FSS. 

At each aperture location, shown in Fig. 3.1, the aperture electric field 
is calculated from the infinite inductive FSS analysis with the plane of the 
FSS positioned tangent to the surface at the aperture center. The necessary 
incident field parameters are the incident magnetic field and the propagation 
vector. The incident quantities are determined by projecting the incident 
vectors found in the previous section into the aperture plane. The procedure is 
repeated for each aperture on the FSS. From the aperture fields, the scattered 
fax fields from the general surface cam be determined. 

3.4. Region III - Forward Scattered Fields 

One method to calculate the transmitted scattered field is by integrating 
the magnetic currents in the spatial domain. This type of field calculation was 
shown by Ko and Mittra in their study of a parabolic cylinder [4]. However, 
since the Fourier transforms of the magnetic currents and aperture electric fields 
are available, there is no need to perform a spatial integration to obtain the far 
fields. The radiated fields can be found conveniently from the transforms of 
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the electric field basis functions. The radiation in a given direction is specified 
by the transmission propagation vector k. Fig. 3.4 illustrates the components 
of the propagation vector in the local aperture plane. Also shown axe the 9' 
and <f > / direction angles. For a given aperture n, the Fourier transform of the 
electric field in the plane of the aperture is 

E^x' j / (^z / » I'y' ) (3.5j 

where (k z > , k y < ) are the transform variables in the local coordinate system 
for the n th aperture. The forward scattered far field expressed in the local 
coordinates is [23] : 

e ~jk 0 R, r / _ _ \ 

E» = j k 0 - - - — a#' ( E x >(k x < , ky > ) cos <f>' + E y > ( k x > , k y > ) sin <f>' ) 

2nR s l \ / ( 3 . 6 ) 

+ a^' cos^' (jS y '(k z ', k y ') cos <f>' — E x >(k x >, ^j,/)sin^'^ j . 

The forward scattered field region is illustrated in Fig. 3.5. The scattered field 
pattern is measured by choosing an angle <f> in the global coordinate system 
and scanning 9 through the desired scan range. The receive probe remains 
pointed at the origin of the global coordinate system throughout the scan. 
The angle between the propagation vector from an aperture n and centerline 
of the horn is 9 p (n). At each scan point, (x 3 ,y 9 , z a ), the received electric field 
is separated into the 9 and <j> global coordinates. The receive probe in the 
experiment was an open ended cylindrical waveguide with an antenna pattern 
modeled by 

E( 9 p ) = E 0 cos n 9 p (3.7) 

Eo is the magnitude of the electric field in the direction of the polarization of 
the probe. Depending on which polarization is desired, the horn is rotated for 
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E* or orientation. The received electric field, E* or E^, is the summation of 
the contributions from each aperture modified by the receive probe’s antenna 
pattern. As a note, if the scan tingles are large and the surface is sufficiently 
curved, there will be a point where some of the apertures will no longer be 
visible. The radiation contribution from these apertures is set equal to zero. 
3.5. - Locally Planar Theory Summary 

The combination of the LPT and the spectral domain FSS analysis is an 
efficient technique for finding the approximate radiation pattern for complex 
FSS geometries. In Chapter 5, the LPT theory will be compared to results 
from an experimental study. The accuracy and applicability of the approach 
will be addressed in terms of incident wavefront and FSS curvature. 
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CHAPTER 4 


FREQUENCY RESPONSE OF THE INDUCTIVE FSS: 
NUMERICAL AND EXPERIMENTAL RESULTS 

4.1 Inductive FSS Design 

The FSS was experimentally studied at the NASA Goddard Space Flight 
Center (GSFC) from 1987 - 1989 as a part of the thesis study effort. The 
results from the studies have proved valuable for comparison with the theories 
developed in Chapters 2 and 3. The experiments included frequency response 
as well as scattered field pattern measurements. In this chapter, the transmis- 
sion coefficient experiments and results for the infinite inductive FSS composed 
of rectangular apertures will be presented. 

In the past, inductive FSSs have been constructed at GSFC by punching 
apertures in an aluminum sheet. In previous years, it was necessary to control 
the punching machine by hand to place each aperture in the proper position. 
Fortunately, a new computer-controlled punching machine was available at the 
shop and this was ideal for fabricating the FSS. The experimental FSS’s overall 
size was chosen to be .6 m square. For the infinite FSS theory to accurately 
model the experimental FSS, a large number of apertures is necessary. Initially, 
S-band (2-4 GHz) was considered. However, at this frequency only about 12 
apertures could fit on the sheet. Moving up in frequency, the range of 9-13 
GHz was chosen because a large number (fa 600) of apertures would fit on the 
FSS sheet and test equipment was readily available in this band. A frequency 
of 11 GHz, the center of the range, was selected as the resonant frequency for 
the experimental FSS. 

The FSS aperture dimensions and spacings were designed by the infinite 
FSS theory presented in Chapter 2 for an incident electric field E x . The FSS 
aperture dimension d, shown in Fig 2.1, was initially chosen to be half a wave- 
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length in size at the resonant frequency. The desired frequency response was 
one which was narrow enough for easy identification of the resonant frequency. 
The aperture size b is the principal means for changing the transmission band- 
width. From the theory, the larger the aperture dimension 6, the larger the 
transmission bandwidth. The FSS was also designed to place the Wood’s 
Anomaly well above the resonant frequency. In this way, a single transmitted 
plane wave could be measured to characterize the frequency response of the 
surface. The Wood’s Anomaly frequency depends primarily on the aperture 
to aperture spacing a. The smaller the value of a, the higher the frequency of 
the Wood’s Anomaly. The choice of a depended primarily on leaving enough 
metal between the apertures for structural integrity. A standard punch size 
of 12.7 x 4.23mm (| x j^inch) was finally chosen for the FSS aperture to 
meet the desired response characteristics. If a non-standard hole size was re- 
quired, each aperture would be nibbled out with a small punch leaving an 
undesirable rough scalloped edge. The resonant frequency was fine tuned by 
varying the dimension c. The dimensions of the apertures and spacings for the 
experimental inductive FSSare shown in Fig. 4.1. 

The FSS thickness needed to be as small as possible to approximate the 
zero thickness theoretical surface. The thinnest aluminum sheet which could be 
punched sucessfully on the punching machine was .5 mm thick. This thickness 
corresponds to where A is the resonant frequency. Thus the experimental 
FSS is a good approximation to the theoretical surface. In later experiments, 
the thin FSS proved its versatility since it could be shaped into a hyperbolic 
cylinder and later be released to spring back to a flat surface. 

For completeness, the initial test setup for measuring both the transmis- 
sion coefficient and scattered fields should be mentioned. The original plan 
was to use the GSFC outdoor antenna range for the FSS tests. Gordon Zee- 
man and Vic Rinker of the University of Maryland designed and constructed a 
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Figure 4.1. Experimental inductive FSS design with ll - GHz resonant fre- 
quency. The dimensions are in millimeters and the FSS thickness is .5 mm. 
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FSS test stand for use with the outdoor range. A number of difficulties such as 
wind, weather, unknown moisture in the wooden stand, and scattering around 
the stand caused the outdoor experiment to be abandoned and the experiment 
was moved inside. 

4.2 Transmission Coefficient Magnitude: Numerical and Experimen- 
tal Results 

Figure 4.2 shows a schematic of the FSS experimental configuration used L 

to measure the magnitude of the transmission coefficient. The FSS is posi- 
tioned at an angle 6 to the incident plane wave direction and surrounded by 
an absorber wall. The orientation of the FSS apertures is given by the co- 
ordinate axes shown. The incident electric field from the transmit horn was 
directed across the narrow dimension of the aperture. The transmit horn was 
placed a distance D — 2.1m away from the FSS and the receive probe was 
placed a distance d = .7 m from the FSS. The receive probe was a waveguide 
to coax adapter surrounded by absorber. The test equipment was configured 
to sweep over the frequency range of interest and the resulting transmission 
coefficient magnitude versus frequency was outputted via a plotter. An initial | 

test was to determine the dynamic range of the measurement setup. By block- 1 

| 

ing the FSS area and comparing the signal level to the case of no blockage, a j 

dynamic range of about 20 dB was achieved. j 

The FSS frequency response was measured ty first calibrating the system 
with the FSS removed. The FSS was placed back in the absorber wall and 
positioned at an angle <5 = 0 for normal plane wave incidence. When the 

frequency sweep was taken, strong reflections occurred between the FSS and | 

1 

the receive probe. The resulting standing wave is seen as an apparent change j 

of the transmission coefficient as the receive probe is moved through a distance j 

range of half a wavelength. Figure 4.3 illustrates the variation of the response 
due to the standing wave. The receive probe distance was moved in 6-mm j 
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RECEIVE PROBE 
AND DETECTOR 



Figure 4.2. Experimental configuration for measuring the magnitude of the 
transmission coefficient. The incident electric field is in the (x,z) plane. The 
FSS orientation is given by the i-axis. The FSS dimensions are found in 
Fig. 4.1. 


43 







10 


10.5 11 11.5 

FREQUENCY GHz 


12 


Figure 4.3. Transmission Coefficient variation due to the standing wave. The 
receive probe was moved in steps of 6 mm about a distance d = .7 m to obtain 
the various curves. The incident electric field is polarized in the x-direction 
with 6 = 0 incidence. The FSS dimensions are found in Fig. 4.1. 
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steps to produce the different curves shown. The standing wave could be 
reduced in magnitude by increasing d because with a fixed receiving antenna 
aperture, the power density at the FSS due to scattering from the receive 
probe is reduced. Shuley has shown the standing wave pertubation can be 
removed from the data by averaging the transmission coefficient as the spacing 
is varied over half a wavelength [24]. This is reasonable if one considers the 
effect of varying d on the response. The standing wave alternates between 
a positive and negative contribution to the receive signal as the distance d is 
varied. Hence the receive power is modulated with respect to distance and can 
be averaged out. 

The solid line in Fig. 4.4 illustrates the averaged experimental transmis- 
sion coefficient obtained from Fig. 4.3. The dotted curve is the calculated 
frequency response for the FSS using 4 full domain basis functions in each 
direction and a large number (101) spectral terms in the inner product sum- 
mations. Convergence of the solution was established with this number of 
spectral terms. The concept of relative convergence has been investigated as 
it applies to the FSS problem. A criterion for the selection of the number of 
spectral terms has been suggested [25]. If the number of basis functions in the 
x and y directions are given by n x and n y respectively, relative convergence 
states M and N should be chosen to satisfy 


n x _ b 
M ~ a 

ny_ = d 

N e’ 

where M and N are defined by 


(4.1a) 

(4.16) 


M = 2m + 1 (4.2a) 

N = 2 n + 1. (4-26). 
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Full domain basis - Large number of spectral terms. 


Figtire 4.4. Experimental and numerical power transmission coefficient as a 
function of frequency. The incident electric field was polarized in the x- 
direction with incident angle 6 = 0. The FSS dimensions are found in Fig. 4.1. 
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The constants a and c axe the unit cell dimensions and b and d are the aper- 
ture dimensions given in Fig. 2.1. The number of basis functions and spectral 
terms axe integer quantities so the equalities in (4.1) and (4.2) can only be ap- 
proximately met. With the experimental FSS geometry and 4 basis functions 
in each direction, the relative convergence criteria states M = 48 and N = 8. 
The numerical result using (4.1) is shown by the dashed line in Fig. 4.4. 

Finally, the subsectional rooftop basis functions were used to solve the 
FSS problem. The basis function set consisted of 3 basis functions in x and 
3 basis functions in y, for a total of 9 basis functions in each direction. The 
resulting matrix equation had dimensions 18 x 18 as compared to 8 x 8 for 
the full-domain case. A laxge number of spectral terms, m = n = 50, was 
chosen for the inner product summations. The numerical results are shown by 
the dot-dashed line in Fig. 4.4. A full discussion of the experimental results 
and numerical calculations has been presented by Webb, Grounds, and the 
author [26]. The experimental resonant frequency was found to be 11.05 GHz. 
The full domain relative convergence solution resonance occurred at about 11.0 
GHz and the laxge number of spectral terms solution resonance occurred at 11.1 
GHz. The subsectional basis function solution also resonated at 11.1 GHz. All of 
the solutions are extremely close to predicting the experimental resonance and 
the accuracy of the experiment could be on the order of the difference between 
theory and experiment. 

Increasing the angle 8 to a non-zero value results in a shifting down of 
the resonant frequency. Figure 4.5 shows the results of the experiment and the 
full- domain, large number of spectral terms, numerical solution. Note how the 
resonance has shifted down to the 10-GHz region. The frequency of the reso- 
nance is quite sensitive to the incident angle 8. In the experiment, the error in 
determining 8 was on the order of 2°. Changing 8 from 30° to 32° produced 
a calculated resonance shift of .175 GHz. Thus, the experimental error could 
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Figure 4.5. Power transmission coefficient for the inductive FSS as a function j 

of frequency. The incident electric field is polarized in the ar-direction with 
8 = 30°. The FSS theory shown is the large number of spectral terms solution. . 

The calculated experimental error is on the order of .175 GHz. The FSS \ 

= •— - 1 

dimensions are found in Fig. 4.1. | 



be a major part of the difference between experiment and theory. 

The experimental FSS discussed represents a very basic filter element. A 
large variety of aperture geometries have been designed to enhance the filter 
characteristics of a FSS. Some of the geometries which have been studied in 
detail are the crossed dipole, Jerusalem cross, and circular aperture [27], [28]. 

In addition, a number of FSSs can be cascaded to produce the RF equivalent 
of a multi-pole filter [29]. For accurate results, the evanescent modes must be 
included in the analysis of the multiple-layered FSSs [30]. 

4.3 Transmission Coefficient Phase: Numerical and Experimental 
Results 

While the magnitude of the transmission coefficient is extremely useful for 
characterizing the overall response of the FSS, the phase response is equally 
as important. For example, if the FSS is used in an antenna system, the FSS 
amplitude and phase characteristics must be included for accurate antenna 
pattern results to be obtained. Just as was the case for the magnitude of the 
transmission coefficient, the phase of the transmitted wave is dependent upon 
both frequency and angle of incidence. 

Figure 4.6 shows the experimental setup for measuring the phase response 
of the inductive FSS. The testing frequency of interest is fed from the signal 
generator to the transmit horn. A sample of the signal is attenuated and mixed 
with a local oscillator signal from the receiver to provide a 45-MHz IF reference 
signal for the receiver. The received signal is picked up by a horn antenna, 
mixed with the local oscillator, and sent to the receiver. At the receiver, the 
phase difference between the received signal and the reference is measured. 
The mixers are located external to the receiver to reduce the large attenuation 
of microwave signals in the coaxial cable. Once the signal is downconverted to 45 
MHz, the coaxial cable loss is small and the receiving antenna can be placed 
at a distant point from the receiver. 
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Figure 4.6. Transmission coefficient phase measurement experimental setup. 
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The phase experiment began by establishing a phase reference between 
the transmit and receive horn with the FSS removed. The FSS was inserted 
between the transmit and receive probes and the phase shift, also called the 
insertion phase, was obtained from a meter on the receiver. The procedure 
was repeated for a number of frequencies to establish the phase response as a 
function of frequency. Figure 4*7 is the measured and calculated phase response 
for the FSS with normal plane wave incidence, 5 = 0. Note at resonance the 
phase shift is equal to zero. 

A large number of FSS problems can be solved with the experimental and 
theoretical results discussed in this chapter given the incident wave is nearly 
planar and the FSS is large and flat. In the next chapter, the experimental 
and numerical results for the LPT will be presented as a way to extend the 
infinite planar theory to general inductive FSS forward scattering problems. 
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Figure 4.7. Transmission coefficient phase experimental and theoretical results 
for a frequency of 11 GHz. The incident electric field is polarized in the x- 
direction and 8 = 0. The FSS dimensions are found in Fig. 4.1. 
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CHAPTER 5 


INDUCTIVE FSS FORWARD SCATTERED FIELD 
EXPERIMENTS AND MODELING 

5.1. Introduction 

Knowledge of the scattered field from a FSS is important for the design 
of FSS subreflectors and radomes. Both the reverse (reflected) and forward 
(transmitted) scattered fields may be necessary for a given FSS problem. In 
this chapter, the forward scattered field is measured and compared with theory 
for a number of FSS and feed/receive probe geometries. The forward scattered 
field was chosen for study because the FSS could be placed in an absorber 
window for good isolation between the feed and receive probes. If the reverse 
scattered field had been chosen, at near normal incidence angles, the receive 
probe would block the illuminating field making it difficult to measure the true 
backscatter response. 

The equipment schematic for measuring the forward scattered field pat- 
tern is shown in Fig. 5.1. The basic testing equipment is identical to the 
equipment used in the Chapter 4 phase measurement study. An antenna po- 
sitioner and chart recorder have been added to complete the test setup. The 
testing frequency, if not mentioned otherwise, is 11 GHz. A photograph of 
the test equipment is shown in Fig. 5.2. In the middle of the left rack is the 
receiver with cables running to the associated equipment. Directly below, on 
the bottom of the left rack is the chart recorder. In the center of the console is 
the antenna positioner and feed horn polarization controls. A number of the 
experiments were performed in a 15 m anechoic chamber when it was available. 
The remainder of experiments were conducted in a large laboratory room with 
microwave absorber blocks situated in appropriate positions to obtain a near 
anechoic environment. 
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Figure 5.1. Schematic of the test equipment used to obtain the scattered field 
measurements. 
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Figure 5.2. Photograph of the Scientific Atlanta antenna pattern measurement 
equipment. 
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5.2. Scattered Field Experiments with Plane Wave Incidence on a 
Planar FSS 

The first scattered field experiment, illustrated in Fig. 5.3, was performed 
in the anechoic chamber. The transmit horn was placed a distance D = 15 m 
from the FSS with polarization shown. The FSS was supported in a large 
absorber wall positioned at an angle 6 from the transmit horn axis. The 
FSS’s orientation is given by the coordinate axes shown. The receive probe, 
embedded in absorber, was placed on the anechoic chamber’s azimuth rotator 
a distance d = .8 m away from the FSS. The receive probe scanned through an 
angle 9 to measure the forward scattered pattern. Figure 5.4 is a photograph 
of the absorber wall in the anechoic chamber positioned at 8 = —30 . The FSS 
has been removed and the receive probe- absorber is visible thru the window 
in the absorber wall. Figure 5.5 is a photograph of the receive probe in the 
forward scattered field area of the experiment. The receive probe is surrounded 
by absorber and the FSS is located in the absorber wall. 

In this experiment, the incident wavefront is nearly planar. Thus each 
aperture’s transmission coefficient is approximately the same. As a first ap- 
proximation, the forward scattered field from the FSS can be calculated by the 
array factor. The array factor pattern is given by a superposition of the fields 
from point sources located at the center of each aperture. The point sources 
are excited with an amplitude and phase corresponding to the incident field at 
the location of the point source. Since the element pattern is broad compared 
to the array factor, the latter dominates the scattered field response. 

The FSS was masked with conducting foil to cover all the apertures except 
for a strip ,3 apertures wide,centered along the y-axis. By choosing a limited 
number of apertures, the sidelobe structure of the array became distinct and 
comparison of theory and experiment was made much easier. Figure 5.6 shows 
the scattered field experimental and numerical results for the case of 3 aper- 
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Figure 5.3. Schematic of the anechoic chamber experiment for plane wave in- 
cidence - planar FSS forward scattered field measurements. The FSS aperture 
orientation is given by Fig. 2.1. 
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Figure 5.4. Absorber wall positioned at 6 = — 30° in the anechoic chamber. 
The FSS has been removed and the receive probe area is visible. 


58 







SCAN ANGLE THETA ( DEGREES ) 


Figure 5.6. Forward scattered field pattern for 8 — 0 plane wave illumination on 
a FSS with 3 apertures in the x-direction and 22 apertures in the redirection. 
The test frequency was 11 GHz, 6 is in the (x,y) plane, and the incident 
electric field is E z , 
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tures in the x-direction and 22 apertures in the y-direction. The FSS was 
positioned at an angle <5 = 0. The experimental scattered field measurement is 
shown by the solid line. The dotted line is the array factor approximation. The 
dashed line is the result using the LPT approach from Chapter 4. The LPT 
takes into account the aperture’s element pattern resulting in a close agreement 
between the LPT results and the experimental data. In all of the scattered 
field experiments, the maximum received power from the experiments and nu- 
merical calculations is normalized to 0 dB. A number of other experiments 
with different aperture masking and orientation angle 6 were explored with 
the array factor approximation method. In general, the array factor predicted 
the scattered field peaks and nulls quite well. The results are presented in 
the paper Finite and Infinite Frequency Selective Surfaces: Experiment and 
Models by Caroglanian et al. [31]. 

In the forward scattering experiment discussed, only a single plane wave 
was incident upon the FSS surface. To explore the forward scattered field 
response of the FSS for plane waves incident at different angles, the experiment 
illustrated in Fig. 5.7 was configured. In this setup, the feed remained the same, 
however, now both the FSS and receive probe rotate as a unit on the azimuth 
positioner with the probe axis perpendicular to the FSS. The forward scattered 
field is measured by rotating the azimuth positioner through an angle 6. The 
phased array model is no longer valid because at each incident angle 0, the 
amplitude and phase response of the FSS is different. Also, since the receive 
probe has been moved up to d = .114 m away from the FSS, the receive probe 
and FSS aperture patterns need to be included in the calculations. The LPT 
method takes into account the receive probe and aperture patterns and is the 
technique studied for the remainder of the chapter. The FSS was reduced in 
size by masking the surface with conducting tape so it was 17 apertures wide 
in the x-direction and 8 apertures wide in the y-direction. The scattered field 


61 


AZIMUTH ROTATOR 



RECEIVE PROBE 


i 

| 

= 

Figure 5.7. Experimental configuration for measuring the scattered field from j 
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the FSS with the incident plane wave sweeping through an angle 6 . Note the FSS « 

and receive probe both rotate on the azimuth positioner. The FSS orientation - 

is given in Fig. 2.1. = 
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Figure 5.8. Forward scattered field as a function of incident angle 9 for a FSS 
with 17 apertures in the x-direction and 8 apertures in the y-direction. The 
operating frequency is 11 GHz. The test setup is shown in Fig. 5.7 with 
D = 15 m and d = .114 m. The FSS dimensions are found in Fig. 4.1 and the 
incident electric field is E x . 
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Figure 5.9. Forward scattered field as a function of incident angle 9 for a FSS 
with 17 apertures in the x-direction and 8 apertures in the y-direction. The 
operating frequency has been changed to 10 GHz. The test setup is shown in 
Fig. 5.7 with D = 15 m and d = ,114 m. The FSS dimensions are found in 
Fig. 4.1 and the incident electric field is E x . 
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pattern measured at 11 GHz is shown by the shaded area in Fig. 5.8. The two 
lines bounding the shaded region are data from the positive and negative 6 
scans. By shading in the difference between the two curves, the non-symmetry 
of the experimental data is revealed. The LPT calculation is shown by the 
dashed line. To check the behavior of the theory with respect to frequency, 
the experiment was repeated at 10 GHz. Figure 5.9 shows the experimental 
and LPT results. Note how the main lobes have flattened and there is a 
transmission peak at 0 = 30°. This peak corresponds to the transmission 
maximum found in the infinite FSS study illustrated in Fig. 4.5 for an incidence 
angle of 30°. 

The experiments in this section have explored the scattered field from the 
planar FSS illuminated by a plane wave. The results in this section were quite 
good because the incident wavefront was nearly planar and the surface was 
flat. However, the LPT can address non-planar incidence as well as curved 
FSSs. In the next two sections, the capabilities and limitations of the LPT 
will be studied in terms of incident wavefront and FSS curvature. 

5.3. Study of the LPT Accuracy as a Function of Incident Field 
Curvature 

The LPT assumes the incident field can be approximated by a number of 
plane waves incident on different regions of the FSS. An important question is 
how does the accuracy of the forward scattered field depend on the incident field 
wavefront? As the incident plane wave becomes non-planar, the infinite FSS 
theory invoked at each subarray becomes a progressively poorer approximation 
to the actual surface response. In this section the limitations of the LPT 
accuracy due to non-planar illumination will be studied by an experimental 
and numerical investigation of the forward scattered fields from a planar FSS 
with spherical wave illumination. The experimental and LPT results will be 
compared as the distance between the feed source and FSS is reduced until a 
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degradation in the numerical analysis is observed. 

Figure 5.10 illustrates the experimental setup for the study. The FSS 
was masked off to be 11 x 8 apertures wide in the x and y-directions, respec- 
tively. The receive probe was positioned at d = .8 m from the FSS. The choice 
of receive probe distance was arbitrary. However, there are two reasons for 
choosing a large value of d. The first reason is to reduce the standing wave 
effect discussed in Chapter 4 for the angles around 0 = 0°. The second rea- 
son is the nulls in the experimental scattered field pattern fill as the receive 
probe is moved closer to the FSS. Thus, the larger d value provides a better 
resolution of the scattered field pattern nulls. This makes it easier to compare 
theory and experiment. The experiment was performed in a large laboratory 
room with absorber blocks placed around the experiment. Figure 5.11 is a 
photograph of the feed source. The feed was a pyramidal horn placed in an 
absorber enclosure. Figure 5.12 shows a view taken from behind the feed horn 
looking at the FSS placed in the absorber wall. The feed was mounted on 
a tripod and was moved to various locations to investigate the effect of dif- 
ferent spherical wavefront curvatures on the experimental and LPT response. 
The masking of the apertures shown in the photo is different from the actual 
masking used in the experiment. Figure 5.13 is a view of the scattered field 
area of the experiment. The receive probe consisted of an open ended circular 
waveguide with a ground plane choke ring surounded by absorber. While not 
obvious from the photo, the rotating positioner’s center of rotation is directly 
below the FSS. Also visible on the absorber wall is a pair of wooden forms 
used to bend the FSS into a hyperbolic cylinder for non-planar testing. 

The experiment began by placing the feed at a distance D = .76 m. Figure 
5.14 shows the experimental and LPT results for D = .76 m. The sidelobe 
peaks and nulls are well predicted for this feed to FSS distance. The feed was 
moved in to D = .46 m and the experiment was repeated. The results, given in 
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Figure 5.10. Experimental configuration for measuring the scattered field from 
the FSS to determine the limitations of the LPT in terms of incident field 
curvature. The distance ID will be varied in the experiment. The receive 
probe distance remains fixed at d = .8 m. 
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Figure 5.11. Photograph of the pyramidal horn transmit antenna. This setup 
was also the receive probe and FSS support for the experiment illustrated in 
Fig. 5.7. 
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Figure 5.12. Photograph of experiment looking from the rear of the transmit 
horn to the FSS positioned in the absorber wall. When the photo was taken 
the FSS was masked for the hyperbolic cylinder tests presented in the next 
section. 
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Figure 5.13. Photograph of the forward scattered field region showing the re- 
ceive probe on the rotating arm. The positioner’s center of rotation is located 
directly below the center of the FSS. Also visible is a wooden template used 
to form the hyperbolic cylinder FSS geometry discussed in the next section. 
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Figure 5.14. Scattered field experimental result for the configuration given in 
Fig. 5.10. The transmit horn is located at D = .76 m. The incident electric 
field is polarized along the x-axis and the frequency is 11 GHz. The FSS 
dimensions axe found in Fig. 4.1. 
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Figure 5.15. Forward scattered field for D = .46 m. At this distance, the LPT 
provides satisfactory prediction of the experimental sidelobe peaks and nulls. 
The incident electric field is polarized along the ar-axis and the frequency is 11 
GHz. The FSS dimensions axe found in Fig. 4.1. 
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Fig. 5.15, show the main lobe and sidelobe structure are still well predicted. 
However, the LPT null at the 30° scan angle is filling. Note at this illumination 
distance the transmission response from 0 - 25° is predicted remarkably well. 
This distance represents the closest the feed could come to the FSS and still 
retain good agreement between the LPT and experiment. Figures 5.16 and 5.17 
show the scattered field results for D = .3 m and D = ,13 m, respectively. At 
D = .3 m the experimental null at 30° is no longer predicted and at D = ,13 m, 
the LPT calculation degrades even further. The failure of the LPT is a gradual 
process and the determination of the distance D is dependent on the degree 
of accuracy desired. 

The loss of accuracy in the LPT is due to the difference in the spherical 
wavefront’s phase distribution over the FSS as compared to the phase of a 
plane wave. To visualize the phase error between the plane and spherical 
wavefronts, consider Fig. 5.18. The transmit horn is positioned at a distance 
D from the FSS represented as the horizontal line. A unit cell of the FSS is 
given by the dimension W» The dotted line P represents the wavefront from 
a plane wave incident on the FSS center. Due to the proximity of the feed, 
there is phase error £ in degrees between the center of the unit cell and the 
edge given by the expression 


£ = 360 (51) 

where \ is the resonant wavelength for the FSS. Also shown in the figure is 
the case of non-normal incidence upon the FSS, illustrated by the subscripted 
variables. By considering the geometry, one can find the phase difference 
between the center and edge of the off-axis unit cell is always less than the 
phase difference for the on-axis case. Figure 5.19 is a graph of the phase error 
between a spherical and plane wave illumination for the experimental unit cell. 
The thin vertical lines represent the distances at which each experiment was 
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Figure 5.16. Forward scattered field for D = .3 m. The 30“ experimental null is 
no longer predicted by the LPT. The incident electric field is polarized along 
the i- axis and the frequency is 11 GHz. The FSS dimensions are found in 
Fig. 4.1. 
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Figure 5.17. Forward scattered field for D — ,13m. Only the first 20 of the 
pattern is predicted. The incident electric field is polarized along the x-axis 
and the frequency is 11 GHz. The FSS dimensions are found in Fig. 4.1. 
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Figure' 5.18. Diagram for calculating the phase error over a unit cell due to a 
spherical source located a distance D from the FSS. 
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Figure 5.19. Phase error over a unit cell due to spherical wave illumination at 
a distance D from the FSS. The solid and dashed curves represent the y and 
^-direction unit cell dimensions, respectively. The resonant wavelength is A. 
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performed. Also shown for each distance is the equivalent distance measured 
in resonant wavelengths A. From the experimental and numerical results, the 
region where the LPT theory is no longer valid occurs at around D = .46 m. 
At closer distances, the phase error rises rapidly resulting in the LPT and 
experimental scattered field differences. So, for the experimental planar FSS, 
the LPT and experimental results are in good agreement until a phase error 
over a unit cell of about 2° or 3° occurs. 

5.4. Study of the LPT Accuracy as a Function of Surface Curvature 

The LPT accuracy is expected to deteriorate with increasing curvature 
since the infinite planar solution for each subarray becomes less accurate. The 
limitation of the LPT due to surface curvature was studied by considering the 
forward scattered field from a set of hyperbolic cylinders of varying curvature. 
The hyperbolic cylinders, illustrated in Fig. 5.20, are independent of y. The 
equation for the hyperbola in the z, z plane, for z as a function of z, is 



Several hyperbolic cylinders identified by the notation H0-H3 were studied. 
The hyperbolic parameters h a and hf, for each curve are given in Table 5.1. 
The hyperbolas were formed by taping the FSS to a wooden form as was 
seen in Fig. 5-12. The FSS was masked off to be 29 x 8 apertures in the 
z and y-directions, respectively. Figure 5.20 also illustrates the hyperbolic 
cylinder experimental configuration. The transmit feed horn was positioned 
at a distance D = 1.75 m to avoid any non-planar wave illumination problems. 
The receive probe distance of d = .8 m remained the same as in the previous 
experiment. 

The first hyperbolic cylinder forward scattered field experiment was for 
the case of zero curvature. The experimental and computed scattered field 
results are shown in Fig. 5.21. As expected, the theory and experiment are in 
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Figure 5.20. Hyperbolic cylinder FSS geometry. H0-H3 represent the different 
FSS curvatures studied. The hyperbolic cylinder is extended parallel to the y- 
axis. The hyperbolic parameters are given in Table 5.1. The FSS orientation 
is given by the x-axis. 


79 



Curve 


HO 


HI 


H2 


H3 


ha 

0 


.217 


.217 


.217 


h b 


.508 


.281 


.178 


Table 5.1. Table of hyperbola parameters for the different hyperbolic cylinders 
studied. The parameter dimensions are in meters. 
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Figure 5.21. Forward scattered field from the hyperbolic cylinder HO. The FSS 
is flat in this case. The incident electric field is polarized along the x-axis and 
the frequency is 11 GHz. The FSS dimensions are found in Fig. 4.1. 
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good agreement. Figure 5.22 is the result for hyperbolic cylinder HI. The main 
lobe was well predicted* however, the sidelobes beyond 30° are not predicted 
by the LPT. The large shaded region is due to the rapid oscillation of the 
sidelobe structure. Figure 5.23 is the case for hyperbolic cylinder H2. The 
main beam is well predicted and the first sidelobe is predicted in position but 
the LPT exhibits a large error in magnitude. Once again after about 30° 
the measured sidelobes sire greater them the LPT calculated values. Finally, 
Fig. 5.24 shows the results from hyperbolic cylinder H3. The curvature in this 
case is becoming large and thus there should be no surprise that indeed the 
theory and experiment differ. Note that even at this curvature, the main beam 
and first null are predicted quite well. 

The hyperbolic cylinder results do not show the gradual decrease in ac- 
curacy as was seen in the non-planar illumination study. Rather, the response 
out to 30°, or the mainbeam in the case of H3, remained good for all curvatures 
studied. At scans greater than about 30°, the LPT was poor in predicting the 
experimental response. 
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Figure 5.22. Forward scattered field from the hyperbolic cylinder HI. The LPT 
and experiment axe in good agreement out to 30 after which the LPT result 
falls below the oscillating experimental results. The incident electric field is 
polarized along the x-axis and the frequency is 11 GHz. The FSS dimensions 
are found in Fig. 4.1. 
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Figure 5.23. Forward scattered field from the hyperbolic cylinder H2. The main 
lobe is predicted well however the experimental sidelobe level is higher than 
the LPT calculation. The incident electric field is polarized along the rr-axis 
and the frequency is 11 GHz. The FSS dimensions are found in Fig. 4.1. 
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Figure 5.24. Forward scattered field from the hyperbolic cylinder H3. The 
main lobe and first sidelobe are predicted well. However the large curvature in 
this case is causing problems with the accuracy of both theory and experiment. 
The incident electric field is polarized along the x-axis and the frequency is 11 
GHz. The FSS dimensions axe found in Fig. 4.1. 
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CHAPTER 6 


CONCLUSIONS 

As a result of the experimental FSS study, the infinite FSS numerical 
solution has been verified for both the amplitude and phase of the transmis- 
sion coefficient. The results show both full domain and subsectional rooftop 
basis functions can accurately predict the frequency response of the FSS. No 
evidence was found to support the use of a relative convergence criterion. 

The forward scattered field experiments revealed for a single plane wave 
incidence on the planar FSS that the transmitted field can be approximated 
by the array factor. If the angle of incidence changes during the scan of 
the planar FSS, the forward scattered field is well predicted by the LPT. 
The limitation of the LPT in terms of incident wavefront and FSS surface 
curvature was determined by comparisons of theory with experimental results. 
For the planar FSS, with spherical wave illumination, the maximum difference 
in phase between the center and edge of the unit cell for good agreement 
between the LPT model and experimental results was found to be about 2 or 
3 degrees. The experiments to determine the effect of surface curvature on the 
LPT accuracy did not show a simple degradation in accuracy as the surface 
was changed to different hyperbolic curvatures. Rather, for all the hyperbolic 
cylinders studied, the LPT predicted the response accurately near the main 
beam followed by large errors in the sidelobes beyond 30°. 

The results from the FSS study can be applied to determine if the LPT can 
be used in the design of a specific FSS systems. If the FSS aperture geometry 
is much different than the experimental FSS used in the experiments, the LPT 
experiments discussed would need to be repeated to determine the appropriate 
guidelines for the particular surface. 
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